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On the various origins of close-in extrasolar planets 
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ABSTRACT 



The extrasolar planets (EPs) so far detected are very different to the planets in 
our own Solar System. Many of them have Jupiter-like masses and close-in orbits (the 
so-called hot planets, HPs), with orbital periods of only a few days. In this paper, 
we present a new statistical analysis of the observed EPs, focusing on the origin of 
the HPs. Among the several HP formation mechanisms proposed so far, the two main 
formation mechanisms are type II migration and scattering. In both cases, planets form 
beyond the so-called snow-line of the protoplanetary disk and then migrate inward 
due to angular momentum and energy exchange with either the protoplanetary disk 
or with companion planets. Although theoretical studies produce a range of observed 
features, no firm correspondence between the observed EPs and models has yet been 
established. In our analysis, by means of principal component analysis and hierarchical 
cluster analysis, we find convincing indications for the existence of two types of HPs, 
whose parameters reflect physical mechanisms of type II migration and scattering. 
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1 TAXONOMY OF HOT PLANETS 

The present EPs database consists of a rather heterogeneous 
sample of planets, showing great variety in all the mea- 
sured quantities. Statistical analysis provides a necessary 
means to find correlations among various physical param- 
eters involved in planetary formation and evolution. Nev- 
ertheless, we caution that statistical analysis may not be 
able to disclose important relationships due to the complex 
and -mostly unknown- interplay of the involved parameters. 
In order to overcome such problems, we performed a global 
statistical analysis of the EPs, using a novel approach. The 
underlying idea is to find groups of similar EPs and to dis- 
tinguish different EP groups on the basis of their diversity. 
In this work, the concept of similarity and diversity among 
EPs is quantified by means of a distance measure in the 
multifold space of physical parameters. This goal has been 
achieved with the aid of principal component analysis and 
hierarchical cluster analysis. In the presen t analysis we fol- 
lowed the procedure used bv lMarchl (|2007f ), and updated in 
iMarchi fc Orto lani (2008). The databasqjused in this paper 
is updated to July 9 th , 2008. We restricted our analysis to 
those EPs having measurements for five input variables, that 
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is: planetary mass (M p ), semimajor axis (a), eccentricity (e), 
stellar mass {Ms) and metallicity ([Fe/H]). Of 308 EPs (in- 
cluding Jupiter), 252 were finally selected for the analysis. 
The purpose of our analysis is to identify planets which are 
similar with respect to the 5-fold space of the input vari- 
ables. As a result, 6 robust EP clusters have been identified. 
Before focusing on HPs, we briefly outline the general na- 
ture of the clusters. 

Cluster CI is characterized by sub- to jovian-like M p , a < 
1 AU, and super-solar [Fe/H]. Cluster C2 has sub-jovian M p , 
a < 1 AU, sub-solar M s and sub-solar [Fe/H]. Both clusters 
CI and C2 also have low mean eccentricity. Cluster C3 is 
the least populated cluster (14 EPs) and probably has no 
strong significance except the fact that it contains many pe- 
culiar EPs that for different reasons have been rejected by 
the other clusters. They are mostly jovian mass planets, hav- 
ing a remarkably super-solar [Fe/H]. Cluster C4 has mostly 
jovian mass planets, a relatively large a, orbiting solar mass 
stars with a widespread [Fe/H] distribution characterized 
by sub-solar values (Jupiter belongs to this cluster). Clus- 
ter C5 has super-jovian mass, a > 1 AU, super-solar M s 
and sub-solar [Fe/H]. The same holds for cluster C6, except 
for its super-solar [Fe/H] and its higher average eccentricity. 
All the input variables have an important role in defining 
the clusters, in particular a and [Fe/H]. Beyond the general 
traits of the solutions, which will be described in more de- 
tail elsewhere, we focus in this paper on the HPs defined as 
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those having a period less than 12 days (see later for further 
comments on this selection). According to the adopted defi- 
nition, 69 HPs are present in the database. The main results 
of the taxonomy is that HPs have been split into two main 
groups, belonging to clusters CI and C2. In addition, there 
are a few outliers placed in other clusters. 

1.1 Two types of HPs 

One of the most interesting features about the HPs, is the 
shape of their mass distribution (see fig. [TJ . It has two peaks. 
These peaks are thought to be real, even if the shape of 
the present distribution is severely affected by observational 
biases that make the discovery of low mass HPs difficult. 
In the literature, the bodies belonging to the lower mass 
peak are referred to as hot neptunes (HNs) while those be- 
longing to the higher mass peak as hot jupiters (HJs). The- 
oretical models also predict a double peaked distribution 
l|Mordasini et al.ll20071 ; llda fc Linll2008T l. where HNs are ex- 
pected to be much more abundant than HJs. 
It is interesting that our taxonomy splits the majority of HPs 
into two different clusters (except for a few outliers). Most 
of the HPs belonging to the peak at log(M p ) ~ are placed 
in cluster CI (fig. HI. HPs of the cluster C2 have a broad and 
natter mass distribution, and it contains most of the HNs. 
Concerning the semimajor axis distribution, similar consid- 
erations hold. HPs of cluster CI are strongly clustered at 
log(a) ~ -1.35 (a ~ 0.045 AU); while those of cluster C2 
have a flatter a distribution (fig. [TJ . The two groups of HPs 
also have clearly distinctive traits in terms of properties of 
stellar metallicity and -to a lesser extent- to stellar mass. 
These variables are important because they account for the 
environment where the planets formed. Figure [2] shows a re- 
markable result: HPs of cluster CI have super-solar [Fe/H], 
while those of C2 have a sub-solar [Fe/H]. Notice that this 
result is also valid for the whole of CI and C2, and not only 
for their HPs. A similar, but less pronounced, result also 
holds for the stellar mass: HPs of CI have mostly super-solar 
M 3 ; while those of C2 have mostly sub-solar M s (fig. O)- 
The HPs of cluster CI and C2 also have other distinctive 
traits. We found some significant intracluster correlations 
which hold for one cluster but not for the other. These are 
the correlations a-e, M p -[Fe/H], M p -M s (see fig.[3J. The 
semimajor axis of the HPs of cluster CI strongly correlate 
with e, while this correlation is absent for cluster C2. On the 
other hand, HPs of cluster C2 exhibit a strong correlation of 
M p -[Fe/H] and M p -M s , while cluster CI does not (fig.EJ. 
These two plots in turn clearly show that HNs of cluster CI 
are very different from the few belonging to cluster C2. The 
latter have remarkable sub-solar [Fe/H] and sub-solar M„. 
Another important point is why some HPs have been placed 
in clusters other than CI and C2. They are: HD73256b, 
HD68988b, HAT-P-7b and HD1 18203b (C3); XO-4b, HAT- 
P-6b and HD162020b (C4); WASP-14b, CoRoT-Exo-3b and 
XO-3b (C5); HAT-P-2b (C6). It is not clear if they are 
real outliers or if they are misplaced by the clustering al- 
gorithm. However, a close look at their properties shows 
that the main characteristic of these HPs is that they have 
high M p and high M s . Among the outliers, the most massive 
(M p ~ 10 Mj) HPs are present, namely HAT-P-2b, WASP- 
14b, CoRoT-Exo-3b, XO-3b and HD162020b (see fig. [TJ left 
panel). Notice that the massive HD41004Bb is placed in- 



stead in C2, due to the very low stellar mass. In the present 
database, the combination of high M p and high M 3 is rather 
unusual, and this suggests that these HPs may be real out- 
liers. Of course, if more HPs having such characteristics are 
discovered, it is possible that they may be grouped into a 
third class of HPs. Finally, it is also important to note that 
previous results are robust to changes in the selection of HP 
periods, at least in the range from 10 d to 30 d. 



2 PHYSICAL INTERPRETATION 



Our present theoretical understanding of the formation 
of HPs is based on two mod els: planetary mi g ration 
within a protoplanetary disk (|Lin fc Papaloizoul ll985T l 
and planet-planet scatte r ing followed by tidal circu l ariza- 
tion jRasio fc Fordl Il996l : IWeidenschilling fc MarzarH Il996l : 
iLin fc Idall 19971 ). Planetary migration has been shown to be 
an efficient mechanism to produce HPs. Migration can occur 
when planetary embryos are still embedded within the disk 
(type I), or when they are large enough to open a gap in 
the disk (type II). According to the present state of the art, 
the migration seems the best candidate for HP formation. 
On the other hand, early works on planet-planet scattering 
showed that the probability for the scattering model t o pro - 
duce HPs was very low. However, iNagasawa et al.l (|2008f ) 
found that the Kozai mechanism enhances the probability 
significantly. Therefore, it is possible that the scattering con- 
tributes to the formation of hot planets as well as type II 
migration, although the latter may be a main channel. 
Type II migration is more efficient for moderate mass plan- 
ets since the planets have to be massive enough to clear 
a gap in the disk, but not too massive to efficiently ex- 
change angular momentum and energy with the disk it- 
self. Previous simulations have shown that the formation 
of HPs via gravitational scattering among planets and the 
subsequent Kozai mechanism combined with tidal dissipa- 
tion is more likely for dynamically active systems of multiple 
planets, typically containing thr ee or more gas giants (e.g. 
M arzari fc Weidenschillinei |2002j) . The formation of many 
iant planets is preferred in high dust surface density disks 
Ida fc Linll2008r ). Dust surface density scales as 10 [Fc/H, M d , 
where Md is the total disk mass. The latter scales, accord- 
ing to theoretica l models, as ~ Mg , where 0.7 < k < 2.2 
(| Vorobvovl I2008T ) . Therefore high stellar masses and high- 
metallicity disks favor the onset of a scattering phase. It 
must also be noted that smaller planets tend to be scat- 
tered inward during the scattering phase and that the tidal 
circularization is more efficient for planets with small mass 
and large physical radius. The detailed final orbital configu- 
ration of the HPs may vary according to several parameters. 
However, there are a number of general features that can be 
outlined. 

In the case of type II migration, the HPs final position is 
close to the location of the disk's inner edge. This is placed 
near the corotation radius (namely the distance where the 
keplerian period matchs the star spin period), where disk 
material accretes onto the stellar poles following the mag- 
netic field lines. The spin period for young stars having 
0.4 Mq < M a < 1.5 Mq may vary con siderably, from 1 d 
up to 20 d or more l|Herbst et al.l 2007). Therefore, taking 
also into account a wide variety of disk parameters, the fi- 
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nal location of HPs formed by type II migration tends to 
be spread, in the range from 0.01 to 0.1 AU. Concerning 
the eccentricity, at the end of the migra tion phase low e 
values are expected. Recent simulations l|Rice et alj 120081 ') 
suggest that, when HPs reach the stellar magnetosphere cav- 
ities, they may further evolve to smaller a. For moderate- 
to- high M p , the e also increases. On the other hand, the 
location of HPs formed by scattering is determined by the 
location where the tidal force is effective. The tidal strength 
depends on several parameters (planetary radius, planetary 
mass, etc) but the final location tends to be in the range 
0.03-0.08 AU, for typical parameter values. The eccentric- 
ity of the inner planet is excited to values close to unity 
by close scattering and via the subsequent Kozai mecha- 
nism. The resultant small pericenter distance enables the 
planet's eccentricity and semi-major axis to be decreased by 
the tidal dissi pation and moderate e ccentricities can remain 
in some cases (|Nagasawa et alj|2008h . We show this in fig. U 
The outcome of the simulations vary according to tide effi- 
ciency. For the simulations shown here , we followed the tidal 
evolut ion of test planets, according to llvanov fc Papaloizoul 
l|2007l ). for a time span of 10 9 — 10 10 yr. A large number 
of planetary radii, planetary masses and stellar masses have 
been considered, choosen at random within the following 
intervals: 0.4 R, T < R p < 2.8 Rj, 0.5 Mj < M p < 3 Mj, 
0.5 Mq < M s < 1.5 Mq. The initial eccentr icity is randomly 
chosen from the distribution obtained by Nagasaw a et al.l 
(2008). Figure [4] (left panel) shows that at the end of the 
scattering and tidal evolution phases the eccentricity and 
semimajor axes are correlated (compared to the observed 
distribution in fig. [3} , while the right panel shows that the 
resulting semimajor axis distribution is peaked. The value 
of the peak depends on the strength of the tide, and for the 
values used here it is peaked at a ~ 0.05 AU, as observed in 
cluster CI (compare with the observed distribution in fig.[l|. 



3 DISCUSSION AND CONCLUSION 

The two main processes of HP formation produce different 
orbital distributions. On the basis of our taxonomic anal- 
ysis, we identify two types of HPs, which may retain the 
footprints of these two different formation processes. In this 
respect, HPs of cluster C2 and CI may have been formed by 
migration and scattering, respectively. This scenario is sup- 
ported by a number of facts. First of all the a distribution 
and the strong a — e relationship for cluster CI, which do not 
hold for cluster C2. Moreover, also the M v distribution of the 
two clusters support this conclusion: higher mass bodies are 
located in cluster CI, indicating more massive protoplane- 
tary disks, which in turn is confirmed by high [Fe/H] and 
high Ms- In this case, the detected HPs are expected to be 
the least massive for each system, since in a multibody scat- 
tering the least massive planets are more effectively pushed 
inward. The orbits of the most massive planets are, how- 
ever, only slightly affected by the scattering phase, therefore 
they tend to stay close to their formation regions and there- 
fore on relatively large a. Although HJs originated via the 
scattering process are expected to be accompanied by outer 
companions, the latter would be beyond the detectable limit 
by present surveys. We find no significant difference between 
CI and C2 about the multiplicity of planets, but we caution 



that this result is affected by low number statistics. Some 
Jupiter-like mass HPs are also present in cluster C2, but in 
this case they may be the most massive bodies formed in 
these systems, given also the moderate-to-low M s and low 
[Fe/H]. Another interesting point is that most of the HNs 
belong to cluster C2. 

Therefore, if we extrapolate directly the percentage of HPs 
belonging to CI and C2 into the efficiency of formation of 
the two processes, we end up with 50% of HPs formed via 
scattering and 30% via migration. The remaining 20% are 
outliers and may be formed either way. These numbers, how- 
ever, have to be taken with caution, since some degree of 
mixing between the two clusters is expected. On the other 
hand, from a theoretical point of view (see discussion in 
previous sections) type II is expected to be more efficient 
in producing close-in EPs. This fact is not in contradiction 
with our findings since it is possible that planets migrat- 
ing inward by type II may stop before becoming HPs. This 
would be the case, for instance, if the gas in the disk dissi- 
pates before the planet reaches the magnetospheric cavity. 
In this respect, it is interesting that many giant EPs exist 
in the range 0.1-3 AU and that their semimajor axis dis- 
tribution is well explain ed by the type II migration model 
l|Schlaufman et al.ll2009l ). 

An alternative scenario is that the two groups of HPs were 
formed by the same process, and the cluster analysis splits 
them on the basis of their diversity in the input variables. In 
this case, the most distinctive variable would be the metal- 
licity. Formation in low metallicity, moderate-to- low M a , en- 
vironments would have produced the HPs of cluster C2. On 
the other hand, high [Fe/H] and moderate-to-high M a would 
have produced the HPs of cluster CI. The combinations of 
these diversities would have also produced the observed dif- 
ferences in the M p and a distribution of the two groups. 
Although this is a possibility, we think the peculiar traits of 
cluster CI and C2 clearly shown the influences of the two 
formation mechanisms. 
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Figure 1. Left panel: Hot planet mass distributions. The current distribution of all HPs shows two peaks (indicated by arrows). 
Right panel: HP semimajor axis distributions. Hot planets of clusters CI and C2 differ considerably in terms of their semimajor axis 
distributions. HPs of C2 have quite a flat semimajor axis distribution, while those of CI seem to be narrowly peaked at log(a) ~ —1.35 
(a ~ 0.045 AU). 




Figure 2. Left panel: HP metallicity distributions. HPs of cluster CI have super-solar [Fe/H], while those of C2 have a sub-solar [Fe/H]. 
Right panel: HP stellar mass distributions. HPs of CI have mostly super-solar M a \ while many HPs in C2 have sub-solar M s . 
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Figure 3. Significant intracluster correlations. Left panel: a — e distribution for cluster CI and C2. The two groups of HPs show a 
remarkable difference: CI shows a very strong, statistically significant, correlation between a and e (2-tailed probability of O.Of %), while 
C2 has no correlation. Moreover, cluster C2 exhibits two statistically significant correlations. They are the M p — [Fe/H] correlation (middle 
panel) and Mp — M s correlation (right panel), respectively with a 2-tailed probability of 1% and 4%. Such correlations do not hold for 
cluster CI. 
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Figure 4. Details of the final orbital configuration obtained in the scattering model. As for the initial conditions, we adopted an uniform 
distribution (in linear scale) of M p , an uniform distribution in planetary radius, and a gaussian M s distribution, peaked at 1 Mq. fnitial 
semimajor axes were choosen at ~ 1.35Mg AU (half of the snow-line distance). The e distribution has been taken according to fig. 9 of 
Nagasawa et al. 2008, where only e > 0.9 is interesting for the tidal evolution (see text for more details). 



